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Abstract 

When imaging studies (e.g. CT) are used to quantify morphological changes in an anatomical structure, it is necessary to 
understand the extent and source of motion which can give imaging artifacts (e.g. blurring or local distortion). The objective 
of this study was to assess the magnitude of esophageal motion due to cardiac motion. We used retrospective 
electrocardiogram-gated contrast-enhanced computed tomography angiography images for this study. The anatomic 
region from the carina to the bottom of the heart was taken at deep-inspiration breath hold with the patients' arms raised 
above their shoulders, in a position similar to that used for radiation therapy. The esophagus was delineated on the diastolic 
phase of cardiac motion, and deformable registration was used to sequentially deform the images in nearest-neighbor 
phases among the 10 cardiac phases, starting from the diastolic phase. Using the 10 deformation fields generated from the 
deformable registration, the magnitude of the extreme displacements was then calculated for each voxel, and the mean 
and maximum displacement was calculated for each computed tomography slice for each patient. The average maximum 
esophageal displacement due to cardiac motion for all patients was 5.8 mm (standard deviation: 1.6 mm, maximum: 
10.0 mm) in the transverse direction. For 21 of 26 patients, the largest esophageal motion was found in the inferior region 
of the heart; for the other patients, esophageal motion was approximately independent of superior-inferior position. The 
esophagus motion was larger at cardiac phases where the electrocardiogram R-wave occurs. In conclusion, the magnitude 
of esophageal motion near the heart due to cardiac motion is similar to that due to other sources of motion, including 
respiratory motion and intra-fraction motion. A larger cardiac motion will result into larger esophagus motion in a cardiac 
cycle. 
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Introduction 

Many groups have investigated the respiratory motion of 
radiation targets (i.e., tumors) and normal tissues in the thoracic 
cavity and the impact of this motion in diagnostic imaging and 
radiation therapy [1,2,3,4,5,6]. Only a few studies have assessed 
the effect of cardiac motion in radiation therapy [1,3,7], and, to 
our knowledge, none have assessed the impact of cardiac motion 
on esophageal motion. The effects from cardiac and respiratory 
motions may not generally be separable as most therapy sessions 
are performed with free-breathing. Also, although respiratory 
motion can be minimized by either breathhold techniques or 
abdominal compression, any motion due to cardiac pulsing wiU 
remain. This may give a blurring of the dose in radiation therapy 
(where treatments are several minutes long); the effect in imaging 
studies will depend on the temporal resolution of the imaging 
study, but could include blurring or local geometric distortion. 
This may be particularly important for longitudinal studies that 
attempt to examine local morphological changes at the voxel level. 

There are several reasons why it may be important to accurately 
image the esophagus - some examples in radiation therapy are 
given here. The esophagus can be an important dose-limiting 
structure in radiation therapy treatment planning, but its dose- 
response characteristics have not been well characterized. That is, 
although it is known that volumes of the esophagus receiving more 
than 40-50 Gy correlate significantly with acute esophagitis [8], 
the impact of the spatial distribution is not yet understood. Some 



studies have shown that circumferential metrics (i.e., the esoph- 
ageal length receiving a certain full-circumference dose) are 
statistically significant [9, 1 0] , but they have not been shown to be 
superior to dose-volume metrics. One likely reason for this is that 
the spatial extent (location) and severity of esophagitis is currendy 
unknown, so easily determining the exact impact of the spatial 
dose distribution on the local tissue is not possible. Several imaging 
techniques could be used to do this on a voxel-by-voxel basis, thus 
allowing the calculation of a true voxel-by-voxel dose-response 
curve. One such technique uses computed tomography (CT) 
images to measure the esophageal swelling that occurs as a result 
of radiation damage [1 1]. Another possible technique is the use of 
positron emission tomography (PET) images. Although this second 
technique has yet to be used to assess esophageal toxicity, it has 
been successfully used to show a relationship between dose and 
standardized uptake values (SUV) for radiation-induced lung 
toxicity [12]. The effects of motion (attributed to respiration or 
some other source) on these techniques is complex, and wiU 
depend on details of the imaging technique used, but could include 
blurring (depending on temporal resolution) or local geometric 
distortions in the esophageal shape. These effects could be 
particularly relevant when trying to assess changes in the 
esophagus (CT or PET, etc.) on a voxel-by-voxel level). Thus, it 
is important to understand the effects of motion prior to using 
imaging (CT or PET) to understand the dose-response of the 
esophagus. Therefore, we initiated this study to assess the impact 
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of cardiac motion on the esophagus as an important step toward 
understanding the different causes of motion uncertainty. 

Materials and Methods 

Patient data 

A total of 26 patients who had received contrast-enhanced CT- 
based coronary angiography with retrospective electrocardio- 
graphic gating had been randomly selected in a previous study 
[7] . The use of these data for the current study was approved by 
the Institutional Review Board at The University of Texas MD 
Anderson Cancer Center (PAl 2-0340), including a specific waiver 
to obtain written consent from the patient for this retrospective 
data review study. The images were taken using a GE LightSpeed 
VCT 64-slice CT scanner (GE Healthcare, Waukesha, WI) and 
deep-inspiration breath hold with patients' arms raised above their 
shoulders, in a position similar to that used for radiation therapy. 
For each patient, 10 CT data sets were reconstructed over a 
cardiac cycle from the carina to the bottom of the heart. These 
images had 1.25-mm slice thickness and 0.35x0.35 mm^ pixel 
size. The images corresponding to the diastolic phases of cardiac 
motion were imported into a Pinnacle'^ treatment planning system 
(Philips Medical Systems, Fitchburg, WI), and the esophagus was 
delineated on a single, consistent phase using standard delineation 
tools. 

Motion analysis 

For each patient the image sets comprised 10 CT data sets 
(To,Tio,..., T90) over a cardiac cycle. (The images are reconstruct- 
ed from multiple cardiac cycles. The phases of the images are over 
a cardiac cycle from 0% to 90% like the ones in 4D-CT.) We used 
an in-house deformable image registration tool to perform the 
motion analysis based on the 10 phases of CT images. An 
accelerated Demons algorithm [13] was apphed to sequentially 
register the 10 CT images and generate 10 displacement fields 
(To->Tio, Tio->T2o, ■■• T9o->To). We validated the deformation 
visually by comparing the deformed images with the original 
image. The 10 displacement fields were concatenated to form the 
motion trajectory for each voxel in a cardiac cycle. The voxel-level 
motion trajectory is illustrated in Figure 1 . The motion magnitudes 
in three directions — left-right (LR), anterior-posterior (AP), and 
superior-inferior (SI) — were computed for each voxel by measur- 
ing the maximum range of the trajectory in aU three directions. 
The total motion magnitude for each voxel was computed as 



where Mlr, Map, and Msi, are the motion magnitudes in the 
LR, AP, and SI directions, respectively. The motion trajectory for 
voxels inside the esophagus contour constitutes the esophageal 
motion. For each patient, and for each CT slice, the mean and 
maximum motions at each direction (LR, AP, and SI) and the total 
motion magnitudes were calculated. Then, for each patient, the 
maximum values of these mean and maximum motions over all 
slices were calculated. Finally, the mean, standard deviation, and 
maximum range of these values were calculated for the overall 
patient population. Note, when calculating the motion in a slice, 
we noticed that the motion in SI direction was not reliable and 
difficult to validate its accuracy. Therefore, we reported the 
motion in the transverse direction without taking into account the 
SI motion. In addition, the esophagus motion of each cardiac 
phase was analyzed, which might be beneficial for determining the 




Figure 1. Illustration of voxel-level motion trajectory. The 

displacement vectors from To->Tio, Tio->T2o, ■•• T9o->To are concat- 
enated together to form this motion trajectory. The motion magnitudes 
in three directions, Mlr, Map, and Msi, were computed by measuring 
the maximum range of the trajectory in all three directions: left-right 
(LR), anterior-posterior (AP), and superior-inferior (SI). 
doi:10.1 371/journal.pone.00891 26.g001 



optimal imaging temporal resolution or acquisition phase in a 
future study. For each patient, the maximum motion of esophagus 
was determined from the 10 displacement fields for each phase, 
and the averaged maximum motion was calculated from the 
overall patient population for comparison among different phases. 

Results 

The maximum esophageal displacement due to cardiac motion 
varied between patients, up to 10 mm in the transverse direction. 
The average maximum esophageal displacement due to cardiac 
motion for all patients was 5.8 mm (standard deviation: 1.6 mm) 
in the transverse direction. The extent of motion is patient specific. 
Table 1 shows the analyzed results for all 26 patients. The mean 
and maximum peak-to-peak motion of the esophagus was a 
function of the superior-inferior position, as shown in the example 
of Figure 2. For 21 of the 26 patients, the largest esophageal 
motion was found in the inferior region. For the remaining 5 
patients, the esophageal peak-to-peak motion was constant (within 
~ 1 mm) throughout the imaged region. When viewing individual 
axial slices, the largest esophageal motion was in the esophageal 
wall closest to the heart or the aorta. Figure 3 showed the averaged 
maximum esophagus motion versus cardiac phase over 26 
patients, including total motion magnitude and motions in LR, 
AP, and SI directions. The motion was found maximum at T80 
with 3.6 mm and minimum at T60 with 2.2 mm. Overall, the 
motion was larger at phases TO, TIO, T80, and T90, where the 
electrocardiogram R-wave occurs. Figure 3 also showed that SI 
direction had slighdy larger motion than LR and AP directions, 
but this might not be true because images used for analysis had 
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Table 1. The analyzed results for 26 patients. 



Mean motion (mm) Maximum motion (mm) 



LR AP Si Axial Total LR AP SI Axial Total 



1 


1.8 


2.9 


5.0 


3.4 


6.1 


5.3 


5.3 


12.7 


6.0 


12.9 


2 


2.0 


1.6 


2.0 


2.5 


3.3 


6.3 


3.4 


6.4 


6.7 


8.4 


3 


1.6 


2.7 


2.9 


3.2 


3.8 


4.0 


6.5 


8.7 


7.0 


9.8 


4 


1.4 


1.3 


2.1 


1.9 


2.8 


3.1 


2.6 


5.8 


3.8 


6.0 


5 


1.5 


1.2 


2.1 


1.9 


2.7 


4.6 


3.8 


4.5 


5.9 


7.1 


6 


1.1 


1.1 


2.5 


1.6 


2.9 


3.0 


3.1 


7.4 


3.7 


7.5 


7 


2.1 


2.2 


2.5 


3.0 


3.5 


4.8 


4.0 


6.1 


5.8 


8.3 


8 


2.0 


1.7 


2.4 


2.5 


3.4 


5.8 


4.5 


4.9 


6.3 


7.7 


9 


1.6 


1.9 


4.0 


2.5 


4.7 


6.2 


5.5 


8.6 


8.2 


10.2 


10 


2.1 


1.9 


3.5 


2.8 


4.5 


6.5 


4.9 


10.4 


7.4 


11.6 


11 


1.4 


1.3 


2.1 


1.8 


2.7 


4.6 


2.9 


6.4 


5.3 


6.9 


12 


1.7 


1.6 


3.4 


2.2 


4.0 


4.1 


3.8 


8.4 


5.0 


8.4 


13 


1.0 


1.3 


1.8 


1.6 


2.4 


3.2 


3.6 


5.0 


4.8 


6.2 


14 


3.1 


2.5 


4.9 


4.0 


6.3 


8.8 


5.6 


9.8 


10.0 


13.7 


15 


1.3 


1.5 


1.9 


2.0 


2.8 


3.8 


4.5 


5.0 


4.8 


5.7 


16 


3.3 


5.3 


4.3 


6.2 


7.6 


6.0 


7.4 


7.1 


8.7 


11.0 


17 


1.5 


1.3 


2.4 


2.0 


3.1 


4.3 


3.6 


8.2 


5.1 


9.4 


18 


2.5 


2.1 


3.9 


3.3 


5.2 


6.6 


4.5 


1 1.3 


6.9 


12.8 


1 9 


1 .8 


1 .6 


2.1 


2.2 


3.1 


3.3 


3.8 


4.0 


4.0 


5.5 


20 


1.4 


1.6 


1.8 


1.9 


2.5 


4.7 


3.4 


3.7 


5.5 


6.0 


21 


1.6 


1.5 


3.7 


2.2 


4.3 


5.1 


3.7 


7.8 


5.4 


9.4 


22 


1.2 


1.5 


1.9 


1.8 


2.6 


2.6 


2.8 


4.0 


3.3 


4.9 


23 


1.7 


1.3 


2.4 


2.2 


3.3 


6.4 


3.8 


49 


6.6 


8.0 


24 


1.6 


1.9 


3.1 


2.5 


40 


3.6 


45 


9.4 


5.5 


9.9 


25 


1.1 


1.1 


1.6 


1.5 


2.2 


3.0 


3.0 


49 


4.0 


5.2 


26 


1.3 


1.4 


2.9 


1.7 


3.4 


3.3 


2.8 


6.1 


3.8 


6.6 


Mean 


1.7 


1.8 


2.8 


2.5 


3.7 


4.7 


41 


7.0 


5.8 


8.4 


SD 


0.6 


0.9 


1.0 


1.0 


1.3 


1.5 


1.2 


2.4 


1.6 


2.5 



For each patient, and for each CT slice, the mean and maximum motions at each direction — left-right (LR), anterior-posterior (AP), and superior-inferior (SI) — and the 
total motion magnitudes were calculated. Then, for each patient, the maximum values of these mean and maximum motions over all slices were calculated and 
presented in this table. The mean and standard deviation (SD) were calculated for all 26 patients as well. Because the motion in SI direction was not reliable and difficult 
to validate, we also reported the motion in the transverse (Axial) direction without taking into account the SI motion. 
doi:l 0.1 371 /journal.pone.00891 26.t001 



1.25 mm resolution in SI dimension but 0.35 mm in LR and AP 
dimensions, thus resulting the motion calculation in SI direction 
was less accurate than LR and AP directions. 

Discussion 

The results of this study showed that parts of the esophagus 
move more than the rest of the esophagus due to cardiac motion. 
To our knowledge, this is the first study to assess the impact of 
cardiac motion on the position of the esophagus. In fact, only a few 
studies have assessed the impact of cardiac motion on radiation 
therapy at any site. Steppenwolde et al. reported that the cardiac 
beat can cause a displacement of 1-4 mm in the position of 
markers in the lung, depending on their location, which ranged 
from 15 mm to 65 mm from the cardiac or aortic wall [1]. 

The esophagus is subject to other sources of motion, including 
respiratory motion and intrafraction motion (e.g., relative shifts in 
position over an interval of several minutes). Pan et al. reported 
that the esophagus clearly moves during breathing, and the 



magnitude depends largely on the individual and on tiie level of 
the esophagus, ranging from 5 mm in the AP and LR directions at 
the thoracic inlet to larger than 10 mm at the gastroesophageal 
junction in these two directions [14]. Several authors have 
evaluated the respiratory motion of esophageal cancers (rather 
than healthy esophagus studied here). Yaremko et al. [15] used 
deformable registration to map the respiratory motion of the gross 
tumor volume for distal esophageal cancer and found that the 
mean tumor motion was 9 mm (95"" percentile: 20 mm) for a total 
motion magnitude. Patel et al. conducted a similar study that 
included several tumors in the proximal two-thirds of the 
esophagus (corresponding closer to the anatomic region we 
assessed in our study), and the total displacements due to 
respiratory motion were less than 4 mm in the transverse direction 
[16]. Cohen et al. evaluated the intrafraction motion of the 
esophagus [4] . They compared CT images taken using a CT-on- 
raUs system before and after radiotherapy and found that 
intrafraction motion of the esophagus in LR and AP directions 
was greater than 5 mm for 1 3 % of cases and greater than 1 0 mm 
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Figure 2. Variation of mean and maximum esophageal motion as a function of CT slice number. 

doi:10.1371/journal.pone.0089126.g002 



for 4% of cases. Thus, for some parts of the esophagus that are in 
the vicinity of tlie heart, the impact of cardiac motion (from our 
study) appears to be comparable to that of respiratory motion 
(compared with the data of Pan et al). Away from the heart, in the 
distal region, respiratory motion likely dominates, although our 
study did not include images away from the heart. 



Conclusion 

In conclusion, cardiac motion results in displacements of 5— 
10 mm for some parts of the esophagus, depending on patient and 
location. This magnitude is similar to that of other sources of 
motion, including respiratory motion and intrafraction motion. A 
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Figure 3. Averaged maximum esophagus motion versus cardiac phase over 26 patients. Error bars for total motion indicate the standard 

deviation over 26 patients. 

doi:1 0.1 371/journal.pone.00891 26.g003 
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larger cardiac motion will result into larger esophagus motion in a 
cardiac cycle. 
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